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As a structural entity of turbulence, hairpin vortices are believed to play a major role 
in developing and sustaining the turbulence process in the near wall region of turbulent 
boundary layers and may be regarded as the simplest conceptual model that can account 
for the essential features of the wall pressure fluctuations. In this work we focus on fully 
developed typical hairpin vortices and estimate the associated surface pressure distributions 
and their corresponding spectra. On the basis of the attached eddy model, we develop a 
representation of the overall surface pressure spectra in terms of the eddy size distribution. 
Instantaneous wavenumber spectra and spatial correlations are readily derivable from this 
representation. The model is validated by comparison of predicted wavenumber spectra 
and cross-correlations with existing empirical models and experimental data. 


I. Introduction 

T he pressure fluctuations induced by turbulent boundary layers have long been of interest as a signifi- 
cant factor in generating vibration, and hence radiated sound, in many practical applications in which 
turbulent flow wets a flexible wall. The need to reduce the noise induced by boundary layers implies a 
requirement for a simple model capable of providing reliable insights into the characteristics of the wall 
pressure fluctuations. 

Turbulent boundary layers exert a shear stress at the wall, and there is a strong connection between 
this shear stress and the structure of turbulent flows in the vicinity of the wall. Much recent empirical and 
numerical research has shown that the region of the flow with high vorticity near the wall appears as a bunch 
of hairpin vortices and longitudinal tubes. Head and Bandyopadhyay * 1 investigated a turbulent boundary 
layer over a wide range of Reynolds number and observed hairpin vortices which are aligned at 45° to the 
wall and convected with a uniform velocity. While these vortices are undergoing a stretching motion, they 
remain straight and group together to form large scale structures, but they do not appear to be interacting 
strongly with each other. In flow visualization studies, Acarlar and Smith 1 ’ confirmed this feature of the 
hairpin vortices. 

We therefore focus attention on fully developed hairpin vortices inclined at an angle of 45° to the wall and 
estimate the associated surface pressure distributions and their corresponding spectra. On the assumption 
that fully developed hairpin vortices are governed by inviscid dynamics, we obtain an exact formulation for 
the stagnation pressure, in terms of a Green function integral along the vortex lines. We then evaluate the 
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surface static pressure by subtracting the dynamic pressure from the results of this formulation applied to 
our vortex geometry. 

On the basis of the attached eddy model, 3 ’ 4 which implies that the form of the wave number spectrum 
can be deduced from the properties of a single eddy, we develop the expressions needed for the surface 
pressure spectra in terms of eddy number-density. This approach draws on flow visualization evidence, 
which indicates that the number of eddies observed in both streamwise and spanwise directions of the flow 
is inversely proportional to their size. The overall wave number spectrum consists of contributions from 
eddies of all sizes, weighted by the number-density. We then validate our model against existing descriptions 
and measurements of turbulent boundary layer pressures. We examine the predicted wave number spectra 
and compare them with those of the Corcos and Efimtsov models 5,6 and also with the results of flight test 
measurements carried out by NASA, Boeing and ANTK Tupolev on the Tu-144LL supersonic aircraft. 


II. Relationship between turbulent structures and wall pressures 


Many analytical approaches linking velocity and 
wall pressure fluctuations in turbulent boundary lay- 
ers have been proposed. One of them, for incom- 
pressible fluctuations, begins with the Poisson equa- 
tion 8,9 

d 2 (ujUj) 


V p = -p- 


dxidxj 


( 1 ) 



Figure 1. Schematic view of the typical hairpin vortex’s 
generation. 


The dominant contribution to the wall pres- 
sure fluctuations in turbulent boundary layers comes 
from Reynolds stresses, and the pressure field in- 
duced by turbulent flows is obtained by solving this 
equation. However, Reynolds stresses cannot be cal- 
culated explicitly until the complete turbulent flow 
is known. The traditional approach consists of de- 
termining models for Reynolds stresses based on em- 
pirical statistics so that the Navier Stokes equation 
can be closed. The characteristics of the wall pres- 
sure fluctuations have been identified through spec- 
tral analysis. Significant contributions to the high 

wave number pressure components are attributed to activity in the near wall region, while contributions to 
the low wave number components are associated with larger structures in the outer layers. 

The long-time evolution of the vorticity field above the wall is probably not deterministic in the usual 
sense and may be regarded as chaotic. Nevertheless, it is clear 10 that vortices are a dominant feature of 
the instantaneous turbulent flow, and furthermore that some events associated with vortices are observed 
to occur repeatedly. Increasingly, the importance of the coherent structures of turbulent boundary layers 
and their contributions to the wall pressure fluctuations is being recognized. These features help in iden- 
tifying characteristics of the surface pressures and their spectra in forming a model for the pressure field 
beneath turbulent boundary layers. Hence a better understanding of the surface pressures induced by typical 
turbulent structures is required. 

Figure 1 presents a schematic view of the typical hairpin vortex’s generation. 10, 11 The initial disturbance 
starts to grow and generate a shape in streamwise direction. Rapidly a vortex head develops, rises from 
the wall, and bends back in the shear flow. However there is no clear explanation for its 45° orientation, 
other than the fact that the maximum mean rate of strain is at 45°. As time increases, the vortex head 
moves farther from the wall while the legs move progressively closer to the surface, and the streamwise 
structures give way to hairpin vortices until the wake region where the structures are predominantly hairpin- 
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like. The original disturbance spreads laterally and produces subsidiary vortices. The development of the 
subsidiary vortices allows the original disturbance to multiply and spread in the spanwise direction through 
an interaction with the background shear flow. The distance between the vortex legs (A) depends on the 
level of the background shear flow. 


III. Numerical formulation 


Here we derive the equations for determining the 
surface pressure field induced by the hairpin vortex, 
which is regarded as two symmetric thin vortex lines 
inclined at 45° to the wall joined by a horizontal 
head (see Figure 2). In order to satisfies the wall 
boundary condition image vortex lines are placed 
symmetrically beneath them. A suitable cut-off dis- 
tance (c) is adopted to avoid the singularity in the 
self-induction integral. 1 "’ Smith and Metzler 13 have 
suggested that the typical vortex has a core diam- 
eter of about 2>Qp/u t and the spacing of the vortex 
pair is 100 u/u T over a range of Reynolds numbers 
in the near- wall regions. Crow 1 " has shown the re- 
lation between the cut-off distance (c) and vortex 
core diameter (cl) is c/d ~ 1/3. The cut-off distance 
is thus taken as proportional to the vortex core di- 
ameter which is also proportional to the distance 
between the two legs. 

We follow the numerical formulation derived by 

Dhanak and Dowling 1 1 for the surface pressure induced due to a developing horse-shoe shaped vortex. By 
using the vector identity, (u ■ V)u = V(|u • u) — u x uj, and introducing the stagnation pressure, p(x, t) = 
! + ^p\u\ 2 , the momentum equation may be written in the form 



Figure 2. Configuration for the hairpin vortex model. 


P 


,du . 9 

p (— — lwxaj = — Vp + pV u 

C/i 


( 2 ) 


The normal component of the momentum equation should satisfy the rigid wall boundary condition: 

= — p(u x u) z -I- pV 2 u z ~ 0 (3) 


dp 

dz 


2=0 


Taking the divergence of each term in Eq. (2) and using the continuity equation (V • u = 0) for an incom- 
pressible flow, a Poisson equation for the stagnation pressure can be obtained: 


V 2 p = — pV ■ (w x u) 

The solution of the Poisson equation can be determined in terms of the Green function 

^ 4-7ij£ — x\ 4-7ij£ — x_ | 4-7T ^ r r ^ 


( 4 ) 


( 5 ) 


where x is the image point of x. and G(£|x) satisfies 

V 2 G = S(i-x) 


( 6 ) 
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with = 0 on £3 = 0 and G(|£ — x\) — > 0 as |£| -> 00 . By using Green’s theorem and these boundary 
conditions, the stagnation pressure is written in terms of a volume integral over the region in which the 
vorticity is nonzero, and a surface integral due to viscous effects. 

p(x) = f V € G(£|s) ■ p(u x u)d 3 £ + [ pGV 2 u 3 d^ 2 (7) 

Jv Js 

If one assumes that the fully developed hairpin vortex is governed by inviscid dynamics, then the equation 
for the stagnation pressure may be simplified further, such that without viscous terms: 

P(x) = p [ V € G(£|x) • (w x u)d 3 £ (8) 

Jv 

IV. Application of the attached eddy model to the surface pressure spectra 

In experimental and theoretical works, 

Perry et al 3,4 showed that an isolated 
hairpin vortex undergoes a stretching 
motion in which the vortex height L in- 
creases uniformly with time and the dis- 
tance A between the two legs of the hair- 
pin vortex at the wall decreases such that 
the product XL remains constant. They 
also showed a geometrically similar ran- 
dom array of hairpin vortices exists in 
the turbulent wall region and all eddy ar- 
rays have the same velocity scale (~ u T ). 

They applied the observation that the 
only velocity scale in the log-law region 
is u T , the friction velocity, to each indi- 
vidual eddy. Thus the typical hairpin vortex, with length scale A, has vorticity ui ~ u T A -1 and circulation 
F ~ u T A. 

Flow visualisation work that they cite indicates that the number of eddies observed in both streamwise 
(x) and spanwise (y) directions decreases like A -1 (see Figure 3). This implies that the number of eddies per 
unit area around a scale A should be four times that around the scale 2A, and this requirement is fulfilled 
by the eddy number density n(A) = N A -3 , with N constant. The physical interpretation of the increase in 
circulation with size is linked to the A -1 distribution of eddies in the streamwise direction, which implies 
combination as they evolve and grow. 

The key point in deriving spectra from the attached eddy model is that the eddies are uncorrelated (on 
average). This implies that the form of a spectrum can be deduced from the properties of a single eddy. 
This approach is applied to velocity auto and cross-spectra by Perry et al, but it can be equally well used 
for pressure spectra. Here, we develop an alternative approach, in terms of eddy number density, to obtain 
expressions needed for the surface pressure spectra. The spectrum is defined by 

/ OO POO 

/ Ri^^e-^e-^d^dp ( 9 ) 

-OO «/— OO 

with the ensemble average of the surface pressure(p), 

v) = p{x, y)p(x + £, y + rj) (10) 

It can be related to the Fourier transform of the pressure field, 

/ oo roo 

/ p{x,y)e~ lk:c e~' lky dxdy (11) 

-00 J —o o 
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via generalized function theory, with the result that 


(2n) 2 5(k X + k X )8(ky + ky)S(k X ,ky ) = P(k x , ky)P(k X , ky) (12) 

We write the pressure field arising from all eddies, 

P(x,y) = '%2Pe(x-x e ,y-y e ,X) (13) 

N e 

where ( x e ,y e ) is the eddy location. The Fourier transform is thus 

P(k x ,k y ) = ^2P e (k x ,k y , \) e - ik **°e- ik *y° (14) 

iV e 


where N e is the number of eddies and P e (k x , k y , A) is the Fourier transform of the pressure field due to an 
individual eddy, and Eq. (12) becomes 

(27T ) 2 d(k X +k X )d(ky+ky)S(k X ,ky) = ^ P g (k X , ky , A) P g {k X , ky A) h ^ ^ ^ ^ ^ ^ (15) 

N e 


Note that the ensemble average has removed all terms involving contributions from different eddies. Then 
the summation is replaced with an integral over eddy location and scale, giving 


(2-7T ) 2 5{k x + k x )5(k v + k y )S(k x , k y ) 

= f x f JZc n(X)P e (k x ,k y ,X)P e (k x ,k y X)e- i ^+ k ^e- i ^y+ k y^dx e dy e dX 


(16) 


Here X m is the scale of the smallest eddy (around 100 v/u T ), and 8 the boundary layer thickness which is the 
scale of the largest eddy. Integration over k x and k y now yields 

r* 

S(k x ,k y ) = J n(X) | P e (k x , k y , A) | 2 dX (17) 

In the same way, the streamwise and spanwise spectrum are respectively expressed as follow: 


nS pOO 

S(k x ,y)= / n(A) | P e (k x ,y, A) | 2 dydX 

•) Am J — OO 

(18) 

pS poo 

S(x,k y ) = / / n(A) | P e (x,k y , A) | 2 dxdX 

•' A m J — OO 

(19) 


V. Non-dimensionalisation and the ‘universal’ spectrum 

The surface pressure due to the eddy is a function of the fluid density p, the vortex strength T, the eddy 
core size d, the length scale A, and the location x, y; that is, 


p e = f(p,T,d,X,x,y) (20) 

We investigated the dependence of the pressure field on the eddy core size and found that it is weak. The 
non-dimensional surface pressure may then be expressed as 

Pe(p,T,X,x,y) = ^p(j) 2 Pe(x/X,y/X) (21) 
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As only the variation in x, y and A is of interest, the dependence on other variables will be suppressed. The 
nondimensional streamwise wavenumber spectrum is 


1 T 2 ~ - - 

O' (k'x : k'y ■ A) — — p—^-P e (k X: k y ') 
where k x = A k x , k y = A k y , x = x/X, y = y/X and 

/ oo roc 

/ P e (x,y)e- ik **e- ik *ydxdy. 

-oo J — OO 


( 22 ) 


(23) 


Substitution into Eq. (17) with the scaling T ~ u T X yields 

_ _ i r s _ _ 

S(k x ,ky) = (-pu 2 T ) 2 J N e X\P e (k x ,k y ) \ 2 d\ (24) 

Finally, one usually wishes to consider the spectrum in non-dimensionalised form. Equation (24) can be 
manipulated to provide the required expression: 

nHvll = f I PeC^ x ,Xk y ) | 2 dX (25) 

{kpu^ys 2 j- Xrn 


where k x = 5k x , k y = 5k y and A = A/A. Firstly we evaluate the ‘universal’ spectrum 


S(Xk x ,Xk y ) =| P e {Xk x ,Xk y ) \ 2 (26) 

for a range of values of the arguments A k x , X k y , and store this in a look-up table. By using a suitable 
quadrature, we then evaluate Eq. (25). Here, whenever the integrand is required for a given value of A, 
we have A k = Xk. hence defining the argument of the ‘universal’ spectrum and therefore the appropriate 
look-up table location. With the same procedure, the non-dimensional streamwise and spanwise wavenumber 
spectrum are respectively expressed as follows: 


1 roc 


S(k x ) 

{\pn 2 T ) 2 6 J Xm )-. 

S (ky ) ^ 


II 


(\pu 2 ys Jx m J- 


I f 


N e | P e (Xk x ,y ) | 2 dydX 


N e | P e (x , A k y ) \ 2 dxdX 


In these cases, the ‘universal’ spectra are respectively defined by 

/ OO 

I P e (Xk x ,y) | 2 dy 

-OO 

and 

/ OO 

| P e (x, Xk y ) 

-OO 


dx. 


(27) 

(28) 


(29) 

(30) 


VI. Wavenumber spectrum model 

On the basis of the Corcos 5 and Efimtsov models, 6 we derive an expression for the wavenumber spectrum, 
evaluate its value and compare with our results. The Corcos model for the spatial correlation between two 
pressures separated by (£, p) in (x, y) is 

S(Z,ri,w) =$(w)e[-^l“«l/ t/c - a ’’l“ ? 'l /votive] ( 31 ) 
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where a x and a y are parameters chosen to yield the best agreement with experiment and <f>(w) is the point 
pressure spectrum. The Efimtsov model considers the dependence of spatial correlation on the boundary 
layer thickness (6), as well as the streamwise (£) and spanwise (rj) spatial separation and it is obtained by 
replacing a x and a y in Eq. (31) by UJ \ lo \ A x and UJ | w | A y respectively. The correlation scales A x and 
A y are 


^■x 

T 

6 



<^2 

+ Sh 2 + ( a 2 /a 3 ) 2 
.- 1 - 1/2 

+ a 2 I for M > 0.9 


+ 


1 - 1/2 
- for M < 0.75, 
r J 


Sh 2 + (o 6 /o 6 ) 


(32) 

(33) 


where M is the Mach number, and Sh = wS/u T . The average values of the empirical constants, 01 ^ 7 , are 
0.1, 0.728, 1.54, 0.77, 548 and 13.5 respectively. We seek to relate the streamwise and spanwise wavenumber 
spectrum to the frequency spectrum: 


S(k x ) = 

— / / S(£, 0 ,w)e- tt -«d£dcj 

J —OO J — OO 

(34) 

S(ky) = 

1 r°° r°° 

— / / S(0, V ,co)e- ik ^d V dco 

J —00 J — OO 

(35) 

The Corcos and Efimtsov models give: 



S(k x ) = 

1 r°° \ x 

7r J _ 00 ^ (w) 1 -)- A 2 (k x — w/U c ) 2 diJ 

(36) 

S(ky ) = 

1 [ $(w) Ay dw 

W- oo 1 } l + (k v A v ) 2 

(37) 


Efimtsov 6 gives the frequency spectrum as a function of the Strouhal number (Sh), which was approximated 
by the flight-test results, 

2 r 

$(w) = ) * 0.01 tt[1 + 0.02(S'/i ) 2/3 ] -1 (38) 

and the relation p rm s = 3.4r 9 then implies an upper frequency limit corresponding to Sh = 3980 in Eq. (36) 
and (37). 


VII. Results and Discussions 

Figure 4 on page 10 illustrates the surface pressure contours of p e in Eq. (21) due to the different eddy 
core sizes, and it shows that the dependence of the pressure field on the eddy core size is sufficiently weak 
that the non-dimensionalised surface pressure can be represented without its contribution. In figure 5 on 
page 11 streamwise and spanwise pressure spectra and their corresponding universal spectra are shown. On 
the basis of the attached eddy model, the streamwise and spanwise wavenumber spectrum for the whole scale 
of eddies are evaluated and compared with the Efimtsov model in Figure 6 on page 12. There have been 
arguments for the spectrum shape, -1 and -5 slope , 9, 15 but its physical explanation is not yet clear. Perry 
et al 4 suggest that there is a range of wavenumbers where the streamwise velocity spectrum 4>(fe) at a given 
height 2 depends on the mean wall shear velocity (u T ), wave number (k), eddy scale (5) and height (z), but 
not viscosity. Furthermore, at the low end of this range, the velocity fluctuations are associated with eddies 
of scale much greater than z, and so 

$(fc) = 6u 2 T gJk6) (39) 
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where gi is a dimensionless function. Conversely, at the higher end of the range one would expect eddies of 
scale 0(z ) to contribute, and that this contribution should be independent of 5. This implies 

$(fc) = zu 2 g 2 {kz) (40) 


If we now postulate an ‘overlap region’ where both of these forms hold, then 

gi(kS) = £ 
g 2 (kz) 5 


(41) 


implying that gi(kd) ~ (fc<5) -1 , g 2 (kz) ~ ( kz ) _1 , and hence that $(fc) ~ k~ l in this region. Perry et al’s 
experiment supports this feature. It is hard to see how these arguments could be applied directly to the wall 
pressures, as there is no height variable z present. However, the other length scale on which the wavenumber 
spectrum clearly depends is \ m , the smallest eddy scale, and one could perhaps argue for dependencies of 
the form 

S(k) = (l/2pu 2 ) 2 Sgi(kd) (42) 

and 

S(k) = {l/2pu 2 ) 2 \ m g 2 {k\ m ) (43) 

at low and high wavenumbers respectively, in which case any overlap region would have -1 slope. In any 
case, the current model clearly predicts a region of -1 slope in both spectra, whereas this is absent from the 
streamwise Efimtsov spectrum. None of the spectra exhibit a distinct region of -5 slope. Figure 7 (a) on page 
12 shows the spatial correlations obtained by inverse Fourier transform compared with the results of flight 
test measurements carried out by NASA, Boeing and ANTK Tupolev on the Tu-144LL supersonic aircraft. 
Here the spatial separations are normalised by different boundary layer thicknesses: the one marked ‘US’ 
is used by Boeing (5 = 0.49 m) and the others marked ‘Russia’ by the Tupolev ( 5 = 0.30 m) for the same 
experimental conditions (distance from aircraft nose = 49.3 m, M = 0.74). The predicted correlations agree 
well with the Tupolev result. Finally, 2-dimensional universal spectrum is provided in Figure 7 (b) on page 
12 . 


VIII. Conclusion 

The present hairpin vortex model shows that it is capable of predicting the properties of the wall pressure 
field, and is therefore a promising candidate for use in exploring features less readily obtainable by other 
methods, eg off-axis correlations. Additional comparisons of streamwise and spanwise wavenumber spectra 
and correlations are presently being made to improve the model. Use of the model for oblique pressure 
correlations will also be investigated. 
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(C) (d) 

Figure 4. Surface pressure contours according to different eddy core sizes: (a) d=10 4 ^-, (b) d=10 3 ^-, (c) 
d=10 2 — and (d) d=10-^. 
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Figure 5. (a) Streamwise, (b) spanwise surface pressure spectra and their universal spectra (c) and (d) 

respectively. 
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Figure 7. (a) Spatial correlations and (b) Stream-spanwise universal spectrum. 
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